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Based on the experiments carried out over the past decade at microgravity conditions, an overview of our 
current knowledge of bubbly and slug flows is presented. The transition from bubble to slug flow, the void 
fraction and the pressure drop are discussed from the data collected in the literature. The transition from 
bubble to slug flow may be predicted by introducing a critical void fraction that depends on the fluid 
properties and the pipe diameter: however, the role of coalescence which controls this transition is not clearly 
understood. The void fraction may be accurately calculated using a drift-flux mo ^ e j^Hfe[! own ^ rom l° ca l 
measurements that the drift of the gas with respect to the mixture is due to the non uniingfeM distribution 
of void fraction. The pressure drop happens to be controlled by the liquid flow for w f hereas for 

slug flow the experimental results show that pressure drops is larger than expect<?|ffl^l^Bj^his study, the 
guidelines for future research in microgravity are given. 

KEYWORDS Two phase flow Microgravity Bubble slug 


INTRODUCTION 

In the past decade several gas-liquid flow experiments have been conducted in micro- 
gravity conditions, using either free fall in drop tower facilities or parabolic flights in 
aircraft. Space applications, such as life support systems and thermal energy produc- 
tion and transport, have raised some classical technical problems that stimulated the 
development of two-phase flow research at microgravity conditions. Among these are 
the prediction of flow pattern, pressure drop, heat transfer and phase fraction in thermo- 
hydraulic systems. Beyond the design of space systems, reduced gravity two-phase 
flows can address some fundamental questions which remain unsolved. On earth, the 
force balance between the two phases is often dominated by the gravitational force. As 
the gravity is suppressed, a new balance comes into play betweeen inertial, viscous and 
interfacial forces so that the mechanics which govern the interactions between phases 
drastically change. Low gravity condition are particularly favourable to the emergence 
of surface tension as dominating force whose role is often small at 1-g conditions. 
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I n this paper we shall focus the discussion on 
related physical issues: 


bubbly and slug flowln^bes and to the 


• What do we know about the transition between bubble and slug flow and what are 
the physical mechanisms to be understood? 

• Is there a drift between phases, which forces are responsible for the drift and what 
are the future directions of research? 

• Does the presence of the other phase modify the basic mechanism of wall friction 
and what are the unresolved questions? 


These questions will be addressed in the following sections. We shall start from the state 
of knowledge that we have reached with the available experimental data and with some 
new data. These data were recently obtained in air- water flow in tubes of 6, 10 and 
19 mm diameters during parabolic flight experiments run in the "Zero-g” Caravelle 
aircraft. The experimental facility is the two-phase flow loop described by Colin and 
Fabre (1995): it is basically the same as the facility described by Colin et al. (1991). 


TRANSITION FROM BUBBLY TO SLUG FLOW 

Among the many experiments carried out in microgravity most are devoted to the 
interface morphology. Indeed the classification of two-phase flow by various patterns, 
although subjective, is easy to accomplish since it only requires a careful observation of 
the flow. In microgravity this observation is made from the pictures taken with high 
speed film or video camera. Several studies with and without phase change have been 
carried out over the past decade. They have been performed with various fluids, in tubes 
of different diameters D and different lengths L during the short periods of aircraft 
parabolic trajectories (between 10 s and 20 s depending upon the plane) or drop tower 
falls (about 2.2 s for the NASA Lewis Research Center's Drop Tower). The fluids used 
were either air and water, or boiling Refrigerant 12 or 1 14. Some of the more recent 
studies are reported in Table I. Different flow pattern have been identified at different 
superficial velocities of liquid j L and gas j G . Let us recall briefly the main results which 
have been put into light. 

At low void fraction, bubbly flow occurs. At high superficial liquid velocity and low 
superficial gas velocity, small bubbles of a few millimetres appear. These bubbles are 
nearly spherical. Their motion is rectilinear with nearly the same velocity. The 
fluctuating motion is weak, contrary to what is observed in 1-g upward flow where 
bubbles of ellipsoidal shape rise with a fluctuating motion. The size of the bubbles is 
mainly controlled by coalescence as suggested by the evolution of the size distribution 
between the inlet and the outlet of the pipe (Colin et al . , 1991) . As the void fraction 
increases, larger bubbles are created. These bubbles move along the tube axis. When 
their size is comparable to the tube diameter, they take an oblate shape controlled by 
the tube size. Following the classification proposed by Dukler et ai (1988) the 
transition from bubbly to slug flow is defined when some bubbles with diameter larger 
than ID or 2D appear. However, this transition is somewhat arbitrary since the larger 
bubbles grow with an increase in void fraction. There is no physical evidence of 
a pattern "bifurcation” (Colin et a /., 1991). 


BUBBLE AND SLUG FLOW AT MICROGRAVITY 


157 



158 


C COLIN et al. 


The large bubbles are the precursors of cylindrical bubbles that appear in slug flow at 
higher void fractions. These long bubbles which have a smooth interface andzt^pheru 
cal shaped nose, are separated by liquid slugs. These slugs contain smaller spherical 
bubbles moving nearly at the same velocity than the cylindrical bubbles. In contrast to 
1-g upward flow these bubbles are not created by gas entrainment at the rear of the 
cylindrical bubbles: they are simply the residue of the initial bubbly flow which is 
injected at the inlet. 

At the gas velocity increases, the liquid slugs decrease in length. When the liquid slugs 
are short enough they collapse. The resulting pattern consisting of lilquid flowing in the 
form of a film at the wall and gas flowing in the centre, is similar to annular flow. But the 
air core sometimes breaks up and frothy slugs, containing many small bubbles, appear. 
This flow pattern which is often called frothy slug-annular flow (Zhao and Rezkallah, 
1993), is a transition between slug flow and annular flow , which occurs at the highest 
superficial velocities of the gas. 

Let us focus on the transition between bubbly and slug flows. Flow pattern maps for 
bubbly and slug flow are shown in Figures 1,2 and 3 for different data sets, using the 
co-ordinates j L and j G . Figure 1 is plotted with the data from Colin and Fabre( 1995) for 
air-water flow in tubes of different size. Figure 2 is plotted with the data from Colin 
et al. (1991) for air-water flow in tube of 40 mm. Figure 3 uses the data from Bousman 
(1995) for air-water flow. 

Several authors have tried to predict the transition between bubble and slug flows, 
the arbitrary criterion being that the transition occurs when some bubbles in the flow 
reach a diameter equal to the tube diameter. The theories may be classified in three 



FIGURE 1 0 g flow patterns in small tubes (Colin and Fabre, 1995): Data i D = 6mm: obubble, • slug 

flows. Data ./ D= 10 mm: □ bubble, ■ slug flows. Data k D= 19 mm: a bubble, ▲ slug flows Transitions: 
bubbly to slug flow: . Equation (2) a t . = 0.2 Equation (2) a f = 0.45. 
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different approaches. In the first one, Eastman et al (1984), Reddy Karri and Mathur 
(1988), proposed to extend to 0-g conditions the correlations used in horizontal flow by 
Weisman et al. (1979) and Taitel and Dukler(1976) or in vertical flow by Weisman and 
Kang (1981) and Taitel et al (1980). For this approach, the magnitude of the gravity 
vector g is replaced by the value of the residual acceleration obtained in microgravity 
experiments. But some of these laws are inconsistent when g becomes very small, In 
a second approach, Lee (1987) and Reinarts (1993) evaluated the different forces acting 
on the mixture for each flow pattern. The transition from bubble to slug flow is assumed 
to occur when the surface tension force becomes greater than the shear stress in the 
liquid phase. In the third approach, Dukler et al (1988), Colin et al (1991) used 
mechanistic models for each flow pattern to predict the transitions. 

The basic mechanism which controls the transition from bubbly to slug flow is the 
coalescence between bubbles. The rate of coalescence depends upon the number of 
collisions which is related to the bubble packing. Duckler et al (1988) pointed out that 
the void fraction of monodispersed bubbles at minimum packing cannot exceed 0.53. 
They found that the transition occurs at a critical void fraction a c of about 0.45 for data 
sets a and b. (Table I). Assuming that in microgravity, water and air are flowing at the 
same average velocity (U L *U G with U L = jj(\ -a)and U G =jJ a) they proposed 
a relation between j G and j L . Colin et al (1991) ignored this assumption and used 
a drift-flux relationship to determine the gas velocity with respect to the mixture: 

~ Co(Jg ^~Jl) (0 

An experimental value of 1.2 was found for C 0 in bubbly and slug flow. With the use of 
Equation 1, the transition may be expressed as: 


Jl — Jg 


1 -Cq«c 

C 0 a c 


( 2 ) 


For air-water flow in tube with inside diameter of 6, 10 and 19 mm (Fig. 1), the 
transition occurs for a critical value of void fraction a c equal to 0.45, as found by other 
authors (Table I). But in 40mm (Fig. 2) and 25.4mm diameter tube (Fig. 3), the 
transition takes place for a critical value a c of about 0.2. Thus for air-water flow, the 
critical void fraction a c depends on the diameter and drops in the range D = 1 9-25 mm. 
What happens if the fluid properties are changed? The experiments of Reinarts (1993) 
with Refrigerant 12 in a tube of 10.5 mm dia., show that the critical void fraction is near 
0.2 (Table I). This proves that the mechanism of coalescence leading to the transition 
from bubbly to slug flow is influenced by the tube diameter and by the fluid properties, 
the refrigerant R 12 having smaller viscosity and surface tension than water. Increasing 
the diameter and decreasing the viscosity seems to promote the formation of the large 
bubbles. This would normally lead to the conclusion that coalescence is promoted 
when the Reynolds number of the liquid Re L =j L D/v L increases. However data 
collected in the range Re L = 7000-20000 in different tube diameters do not support this 
idea. For same Reynolds numbers range, the transition occurs for a c %0.45 in small 
tube whereas it occurs for a c % 0.2 in 40 mm and 25.4 mm diameter tube, suggesting that 
the transition depends of neither j L nor j G . If we accept this conclusion, the value of a c 
would depend on a dimensionless number which does not involve the superficial 
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velocities. A good candidate is built with the fluid properties, i.e. the liquid viscosity v ; , 
the liquid density p L and the surface tension a: 


N 


D “ 


oD 

pi>’l 


(3) 


N d is also equal to Re\jWe , where We is a Weber number defined as We = p r t j\Djrs. It is 
thus suggested that the transitional void fraction a c be a function of N n : y. c —f (N D ). 

This transition was checked against the existing data in Table I. However, as the void 
fraction is missing in some of them, a was replaced by the volumetric quality of the gas, 
defined as: 


* = (4) 

Jl +./g 

In bubbly flow, Equation 1 holds so that it is always possible to express the transitional 
void fraction a c versus the transitional quality x f : 

x c = C 0 x c . (5) 

Thus we may look for the existence of a relationship of the form x c = g(N D ). 

The range of variation of x c at the transition is plotted in Figure 4, for different sets of 
data listed in Table I. Note that each data set is represented by a vertical segment, since 
for each set, N D is constant. This segment delineates the range within which a transi- 
tional regime has been found at different velocities: it may be interpreted either as the 
range of uncertainty of the observations or as the result of a small variation of x c with 
the mixture velocity. Excepted for the data of Huckerby et «/., the values of x ( . fall 



FIGURE 4 0-g dimensionless flow pattern map: (ZZZZZ3 transition. 
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mostly within the range 0.44-0.54, corresponding to a c between 0.37 and 0.45, for N D 
less than 1.3 x 10 b . For higher values of N D , x c ranges between 0.22 and 0.27, 

Despite the uncertainty in determining the transition, two different regimes are put 
into light. For values of N D below 1.5 x 10 6 , coalescence seems hindered and the 
transition from bubbly to slug flow appears sensitive to the bubble packing at about 
0.45 of void fraction, whereas for values of N D greater than 1.7 x 10 6 , the mechanism of 
coalescence is promoted and leads to a rapid transition to slug flow at a c equal to 0.2. 
An illustration is given in Figure 5 and 6 with the pictures of bubbly flow at both inlet 
and outlet of the tube for the same flow conditions and different tube diameters. For the 
smaller tube (Fig. 5) it appears clearly that the flow pattern does not change drastically 
between inlet and outlet. For the larger tube (Fig. 6) big bubbles are observed at the 
outlet indicating that the collisions between bubbles were very effective in promoting 
coalescence. 

The available results for the transition from bubbly to slug flow may be summarized 
as follows: 


aD 


,< 1.5 x 10 6 =>« «0.45, 
Pi A 

x c % 0.54 

a ° 7 > 1.7 x 10 6 => a «0.20, 

x c % 0.24 

PlK 



This relationship has the advantage to take into account the fluid properties in the 
evolution from bubbly to slug flow, — evolution being here more appropriate than 




FIGURE 6 Flow visualization in 40 mm dia. tube: / L — 0.41 m/s, j 0 =0.1 1 m/s. Left: inlet test section. 
Right: outlet test section (80 D downstream). 
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transition which evokes a bifurcation in flow pattern — . Equation 6 raises an interest- 
ing question: is there a transition between the hindering-coalescence regime and the 
promoting-coalescence regime? The existence of such a transition has to be confirmed 
by additional experiments in the range of N D between 1.5 x 10 6 and 1.7 x 10 6 . It also 
has to be explained based on physical arguments. 

Due to the lack of experimental results these conclusions are still questionable. They 
suggest however the line to follow for future experiments, concerning the transition. 
With small tube, experiments without phase change, the gap will be difficult to fill. 
Indeed the only way would be to reduce the viscosity, but that would use fluids as 
alcohols or solvent, which would conflict with aircraft safety restrictions. In contrast 
experiments in large tubes will be highly valuable since the viscosity can be easily 
increased by addition of glycerine in water and the surface tension can be reduced with 
the use of mineral oil. 


PHASE FRACTION 

Several authors have published data concerning either void fraction or average 
gas velocity. Different methods have been used to determine the cross-sectional 
averaged void fraction a. In the studies of Colin et al. (1991) and Colin and 
Fabre"(1995) conductance probes made of 2 stainless steel flush-mounted electrodes 
inside the tube wall were used. Bousman and Dukler (1993) have used conductance 
probes with two parallel wire in order to determine the void fraction in the bubbly 
region of bubbly and slug flow and the film thickness in the separated flow region 
of slug and annular flow. The cross-sectional averaged gas velocity U G was determined 
from flow pictures (Colin and Fabre, 1995) or by cross correlating the signals from 
two pairs of wire positioned along the tube axis (Bousman and Dukler, 1993). The 
measured values of gas velocities U G were compared to those calculated from the 
experimental void fraction U G — j G j a (Colin and Fabre, 1995). In most of the cases 
a good agreement is found for bubbly flow, in contrast to slug flow for which some 
discrepancies have been observed in small tubes. This is partly due to th£ difficulty to 
calibrating the conductance probes in small tubes, since the apparent electrical 
conductance is strongly sensitive to the local void distribution inside the tube, i.e. to the 
flow pattern. 

The gas velocity U G is plotted versus the mixture velocity) =j G +) L , in Figure 7 for 
bubbly flow and in Figure 8 for slug flow. In Bousman and Dukler's results U G was 
determined from the cross correlation between the signals of two void fraction probes 
located at certain distance from each other whereas in our present results, U G was 
obtained from an image processing method. These results lead to the following 
conclusions: as C 0 > 1 there is an average drift of the gas with respect to the mixture 
equal to U G -j = (C 0 - \)j; regarding the drift, bubbly and slug flow can be hardly 
distinguished. This drift is reasonably well predicted by Equation 1, as already 
proposed by Colin et al (1991): 

U (} = C 0 j with C 0 %1.1-1.3. 

Where does the drift come from in the absence of gravity? 


( 7 ) 
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The local velocity difference is expected to be weak in the absence of body force. This 
is confirmed by the experiments of Kamp et al. (1993) whose results are plotted in 
Figure 9. In this figure the local liquid velociy u L , measured with a hot film anemometer, 
is compared to the local gas velocity u c obtained from a double optical probe. The local 
velocity difference between phase is small: it remains within the range of uncertainty of 
the measurements. These results provide a strong argument that no local slip exists in 
microgravity, thus u G % u, % u. 

Nevertheless, even in the absence of local slip, an average drift can exist if the gas is 
mainly concentrated in the flow regions where the velocity is the greatest. Let e be the 
local void fraction and ( — } = A~ l $ A — dA the averaging operator over the cross 
section A of the tube. Then C 0 may be expressed as: 


Co — 


<eu> 

<*j ’ 


( 8 ) 


C 0 characterizes the relation between the distributions of local velocity and void 
fraction. Thus, the value of C 0 gives some information on the local distribution of void 
fraction. It is greater than unity if the radial distribution of void fraction is at its 
maximum at the tube axis. Kamp et al. (1993) have determined the distribution of local 
void fraction by using optical probes: they found indeed a maximum of the local void 
fraction at the tube axis (Fig. 9). Following Zuber and Findlay (1965) it is possible to 
determine the value of the coefficient C 0 from Equation 8. Assuming that the local 



FIGURE 9 Void fraction profile and velocity profiles of liquid gas in 40mm dia. tube: J L = 0.77 m/s, 
j G = 0.044 m/s. o velocity of liquid, □ velocity of gas, ■ local void fraction, Equation (10) with m = 9; 

Equation (9) with n = 2. 
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velocity and void fraction distributions have the form: 



the coefficient C 0 becomes: 


“="»[' “(s) ] 

m + n + 4 
0 m + « + 2 


( 10 ) 

( 11 ) 


For the test corresponding to j L = 0.77 m/s and j G = 0.044 m/s, the values of m and « are 
found to be equal to 9 and 2 respectively from using a least square fit. The resulting 
value for C 0 is 1.15, in good agreement with the value calculated from the experimental 
void fraction given by the conductance probe. 

The effect of spatial distribution does explain that there exists an average drift in bubbly 
flow even if there is no local slip. Therefore the key issue for bubbly flow is to understand 
the lateral force balance which is responsible for the concentration of void near the axis. 
The answer is not clearly known at the present time and the local two-fluid models fail 
to predict the correct distribution. An effort has to made on this simple and fundamen- 
tal case. Some elementary experiments like the motion of a single bubble in pipe at 
microgravity conditions would be helpful to understand the phenomenon: it could be 
a physical or numerical experiment. On the other hand the experimental analysis of the 
local distributions of velocity and void are also great help. The first experiments carried 
out by Kamp et al. on the determination of the local structure of the flow, show 
however some limitations in performing such experiments in parabolic flights because 
of the short duration and the poor quality of the microgravity conditions in parabolic 
flights. This pleads strongly for an in-orbit experiment which could be decisive. 

As previously mentioned C 0 has for slug flow the same trend than for bubbly flow 
(Fig. 8). This was already pointed out by Colin et al. (1991), based on their experiments 
in 40 mm diameter tube; this is also confirmed by the present experiments in smaller 
tubes and by those of Bousman and Dukler (1993). In slug flow the gas is mostly 
contained in the large cylindrical bubbles. It is thus acceptable to assume that U G % V , 
where V is the velocity of the cylindrical bubbles. The result on C 0 shows that the 
Nicklin et al. (1962) relationship for turbulent flow' V= 1.2 j + 035(gD) 1 2 , is valid also 
in the limit*/ = 0. From the inviscid flow theory of Collins et al. (1978), oneshoul be able 
to predict the value of C 0 for both laminar and turbulent flows of the continuous phase 
in the limit g = 0: 


for laminar flow 


log Re L -I- 0.089 
log Re L — 0.74 


for turbulent flow 


The turbulent flow relationship gives a value decreasing from 1 .25 to 1.2 for Re L in the 
range 10,000-100,000. However this variation is within the range of uncertainty of the 
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measurements. On the other hand at smaller Reynolds number, typically < 5,000, C 0 is 
likely to increase up to reach the theoretical value of 2.27 when the flow becomes 
laminar. In \-g flow this behavior was clearly shown by Frechou (see Fabre and Line, 
1992). Some additional experiments are thus necessary to determine the behavior of C 0 
for Re, within the range 500-5,000. Indeed in this range, the transition -which must 
exist-is highly indicative of a modification of flow regime upstream the nose of the 
cylindrical bubble. 


WALL FRICTION 


In microgravity the friction at the wall is directly related to the pressure drop dP/dx. 
Indeed the mean shear stress t„. may be expressed by 


DdP 

w ~4~dx' 


(13) 


Measuring the pressure drop gives access to the wall shear stress. This holds for 1-g 
horizontal flow as well. Most of the results on pressure drop reported in the literature 
concerns two-phase flow at high void fraction in annular or slug pattern (Chen et al . , 
1991; Bousman and McQuillen, 1994; Zhao and Rezkallah, 1994). Only scarce results 
have been published for bubbly or slug flow. 

It has already been pointed out by Chen et al. (1991) that pressure drops are 
generally greater in microgravity than in 1 -g horizontal flow. In order to illustrate this 
peculiarity, the values of t w , which are calculated from different data sets ( Bousman and 
Dukler, 1993; Colin and Fabre, 1995) have been replotted in Figures 10 and 1 1 for two 
different tube diameters. In Colin and Fabre experiments, the pressure drop is 
measured from two different transducers and only the results which agree within 10% 
are retained. The difference between 0-g and 1-g results is shown to be smaller for the 
smaller tube. This difference also decreases when the gas velocity increases. It is 
expected that, when inertial force is large enough compared to gravity, there is no 
reason to distinguish 0-g flow and 1-g horizontal flow: this is probably true when 
p G U G /ApgD» 1. i.e. in annular flow in small tubes. 

What controls the wall friction in bubbly flow? One method to answer to this 
question goes through the analysis of the friction factor, generally defined as 
f=zJ(\/2)pU 2 . This definition raises the following question: what density and 
velocity scales must be chosen? We have the choice between the liquid quantities p , , L , 
or the mixture quantities p M = ot L p L + x G p G ,j or a combination of both. It is possible for 
example to introduce the two following definitions: 




- 2 P,ul 


(15) 


The fundamental issue is not only a matter of definition, it is intimately related to the 
physics of momentum transfer at the wall. A piece of the answer can be given from 
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FIGURE 12 Wall friction factors f M in bubbly flow: 0-g flow o D = 6mm, □ D = 10mm, a D = 19 mm, 
■ D = 12.7 mm, o D = 40 mm, Blasius relationship 



K‘l 


FIGURE 13 Wall friction factors f L in bubbly flow: 0-g flow o D = 6mm, □ D — 10mm, a D = 19mm, 
■ D = 12.7 mm, o D — 40 mm single-phase flow: + D = 6 mm - D — 10mm —D = 19 mm, Blasius, 

Foiseuille relationship. 
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intuitive arguments: as the void is mainly concentrated at the tube axis, the momentum 
transfer at the wall must be mainly controlled by the liquid motion. It is thus 
anticipated that the friction factor defined by Equation 14 be given under the form 
f L ~ F(Re L , 0 L ), where/) = F(Re L , 0) is the single phase flow relationship. An alternative 
is to express the friction factor defined by Equation 15 as f M = F(Re M , a). The 
experimental results are plotted in Figures 12 and 13 with the two choices of dimen- 
sionless variables: the results are less scattered in Figure 13, justifying that the liquid 
motion controls the shear stress at the wall. In addition, the experimental points are not 
too far from the single phase Blasius relationship: 

f L = 0.019 Re h ° 2 \ (16) 

This relationship may be used with confidence to predict the pressure drop in the range 
Re z = 20,000-80,000. Indeed for bubbly flow in 40 mm diameter tube, it has been 
shown by Kamp et ai (1993) at Re L = 40,000, that the near wall velocity profiles follow 
the same logarithmic law as far single-phase flow: the presence of the bubbles does not 
significantly modify the local structure of the flow. In contrast, for Reynolds number 
less than 20,000, the experimental friction factor is greater than predicted by Equa- 
tion (16). This tendency may be explained as follows: as the Reynolds number de- 
creases, the thickness of the viscous layer increases but the presence of large bubbles 
may affect this layer near the tube wall. This is due to an additional production of 
turbulence induced by the distortion of the liquid flow by the bubble motion: this extra 
turbulence has indeed been put into light in bubbly flow at microgravity conditions by 
Kamp et ai (1993). 

Another interesting result which appears in Figure 13 concerns the transition 
between laminar and turbulent regime. Single phase flow experiments have been 
carried out in the same loop. The corresponding friction factor has been added in 
Figure 1 3. For experiments in 6 and 10 mm tube, the Reynolds number ranges between 
1 ,000 and 40,000. In the 6 mm diameter tube, when the single-phase flow is laminar, the 
wall friction factor follows the theoretical Poiseuille relationship/; — 16/Rc, . The tran- 
sition from laminar to turbulent flow is observed at Re % 8,000 in single-phase flow in 
the 10 mm diameter tube, the jump of f L near Re { = 8,000 being the signature of this 
transition. Clearly two-phase bubbly flow does not display the same behavior at least 
for Re L > 4,000. There are some interesting and yet unresolved questions related to the 
wall friction behavior at smaller Reynolds number. Does a transition exist between 
laminar and turbulent flow and does the wall friction decrease at smaller Reynolds 
numbers? This motivate the needs for careful experiments at low Reynolds number. 

Let us now discuss the case of slug flow from the experimental results plotted in Figure 
14. In comparison to bubbly flow the same trends are observed, although the experimen- 
tal data are much more scattered. There is a fundamental difference between bubbly 
and slug flow: the intermittence. In a reference that moves with the long cylindrical 
bubbles, the flow appears almost frozen in microgravity. Indeed each long bubble is 
separated by bubbly liquid slugs in which the small bubles move at about the same 
velocity as the long ones. By using the unit cell concept proposed by Wallis (1969), the 
pressure drop may be written under the form (see for example Dukler and Fabre, 1 992): 

dP _ Pi^LwB^LwB T t GwB^Gwe) + 0 ~ P)( T LwS$ Lws) 

dx A 


( 17 ) 
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FIGURE 14 Wall friction factors f L in slug flow: 0-g flow o D = 6mm, □ D = 10mm, A D= 19mm, 
u D = 12.7 mm, o D = 40 mm single-phase flow: + D = 6 mm, D = 10mm, D — 19mm, Blasius, 

Poiseuille relationship. 


where the subscripts B and S stand for the long bubble and slug regions respectively, 
A is the crossed sectional area, S is the wetted perimeter and /? the ratio of the long 
bubble length to the total cell length. As the gas does not wet or contact the wall in 
microgravity, the mean shear stress defined by Equation 13 becomes: 

K, = Pt wS + (\-P)z wD . (18) 

For Equations 1 7 and 1 8 to be written under this form one has to accept that the flow is 
fully developed in both regions. In the absence of gravity, there is no driving force to 
move the liquid film around the long cylindrical bubbles: this is verified from video 
sequences, the liquid velocity being estimated from the motion of tiny bubbles in the 
liquid film around the long bubbles. Thus the shear stress in the long bubbles is 
probably weak in contrast to this in the liquid slugs. Consequently it may be assumed 
that t„, «(1 — P) t hd- With the fully developed flow assumption, the velocity of the 
mixture in the liquid slugs is equal to j. If we accept that the friction factor in bubbly 
slugs is given by Equation 16, as shown by the foregoing results, then the shear stress in 
slug flow must be equal to the shear stress in the corresponding bubbly flow, times 
1 - fi. This assumption would lead to the result that the shear stress in slug flow must be 
smaller than the shear stress in the corresponding bubbly flow. This is not supported by 
the experimental results which have shown the opposite. To arrive to this paradox we 
have used two assumptions: the flow in each region is fully developed and the friction in 
the liquid film is weak enough to be disregarded. It is likely that the second assumption 
is correct. This is not the case for the first assumption since the slugs are short and the 
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velocity profile can hardly be considered as fully developed in the moving frame. This 
assumption has to be reconsidered if we want to correctly predict the friction at the wall 
in slug flow. 

CONCLUSION 

During the past decade, a considerable effort has been done to understand two-phase 
flow at microgravity conditions. The studies have generally concerned flow pattern 
transitions and at a least extend void distribution and pressure drop. Focusing our 
attention on bubbly and slug flows, we attempted to present the state of our current 
knowledge from the current available data. Some conclusions may be drawn for the 
future research. 

Concerning the transition from bubble to slug flow, we suggest the existence of 
a possible transition for some value of oDj'pj v ] ; between 1.5 x 10 6 and 1.7 x 10 6 . Future 
experiments have to be carried out to confirm this suggestion. Moreover this transition 
poses the delicate problem of predicting bubble coalescence in flow conditions, which 
must be analysed from this perspective. 

The void fraction may be predicted with an acceptable accuracy from Equation 1. 
This equation proves that a mean drift between gas and liquid exists. For bubbly flow, 
the origin of the drift must be found in the concentration of void at the tube axis. At the 
present time the force balance which controls the bubble migration towards the axis 
remains unclear and specific experiments are needed. For slug flow the motion of long 
bubbles is well explained by the theory in the limit g = 0. However some experiments 
are needed at lower Reynolds numbers, say between 500 and 5,000 to clarify this motion. 

At microgravity conditions the pressure drop is the direct signature of the mean 
shear stress at the wall. In bubbly flow, the single-phase flow relation of Blasius is 
acceptable although there is a tendency to underestimate the friction when Re L < 
40,000. The single-phase relationship ignores the presence of bubbles which needs to be 
taken into account. Some new’ insight into the local mechanism of bubble turbulence 
interaction is needed. For slug flow, the friction is not predicted with a sufficient 
accuracy from the unit cell concept: the fully developed flow assumption has to be 
reconsidered. Last but not least, the evolution of the friction factor with the Reynolds 
number may indicate whether a transition from laminar to turbulent flow exists in 
two-phase flow. In this view, the analysis of the velocity fluctuations may be useless 
because the bubbles induce small velocity fluctuations by their motion relative to the 
wall that may be difficult to discriminate from turbulence fluctuations. 
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